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Abstract 



Following the basic prescriptions of the Tsallis' nonextensive thermodynamics, we 
study the relativistic nonextensive thermodynamics and the equation of state for a 
perfect gas at the equilibrium. The obtained results are used to study the relativistic 
nuclear equation of state in the hadronic and in the quark-gluon plasma phase. We 
show that small deviations from the standard extensive statistics imply remarkable 
effects into the shape of the equation of state. 
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1 Introduction 



It has been shown that the nonextensive generalization of the Boltzmann- 
Gibbs thermostatistics, proposed by Tsallis, can be very relevant in many 
physical applications where long-range interactions, long-range microscopic 
memories and/or fractal space-time constraints are present [1]. In particu- 
lar, recently, many authors outline the possible connection to the nonexten- 
sive statistical framework with nuclear and high energy physical applications 
[2,3,4,5,6,7]. The aim of this work is to generalize the basic concepts of the 
nonextensive statistical mechanics to the relativistic regime and to investigate, 
through the obtained relativistic thermodynamic relations, the relevance of 
nonextensive statistical effects on the hadronic and on the quark-gluon plasma 
(QGP) equation of state (EOS). As we will see, small deviations from the ex- 
tensive thermostatistics produce a significant modification into the shape of 
the hadronic and QGP equation of state with important consequence on the 
deconfined phase transition and on several nuclear properties. 
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2 Relativistic nonextensive thermodynamics 



In this section we present the basic macroscopic thermodynamics variables in 
the language of the nonextensive relativistic kinetic theory. Let us start by 
introducing the particle four-flow in the phase space as [8] 

n"(x) = t f^riTffrp), (l) 



ZnJ 



and the energy-momentum flow as 

T'*( X ) = l r ri2.jfjrF( XlP ) 1 ( 2 ) 

Zj n J p 



where we have set K — c — 1, x = x M = (i, x), p = p^ = (p°,p), p° = 
VP 2 + m' 2 is the relativistic energy and f(x,p) is the particle distribution 
function. The four- vector iV M = (n,j) contains the probability density n = 
n(x) (which is normalized to unity) and the probability flow j = j(x). The 
energy-momentum tensor contains the normalized g-mean expectation value 
of the energy density, as well as the energy flow, the momentum and the 
momentum flow per particle. Its expression follows directly from the definition 
of the mean g-expectation value in nonextensive statistics [1]; for this reason 
it is given in terms of f q (x,p). 

Furthermore, in the framework of the nonextensive thermostatistics, it appears 
natural to generalize the nonextensive four- flow entropy Sg(x) as follows 

W = -K I ^P^ <? (x,p)K/(a;,p) - 1] . (3) 



It is possible to show [8] that such an entropy, together a generalized relativistic 
Boltzmann equation, satisfies the relativistic local if-theorem and implies the 
following Tsallis-like equilibrium probability distribution 



feqip) = tt 



1 - (1 - q) 



k B T 



1/(1-9) 



(4) 



where is the hydrodynamic four-velocity [9] and f eq depends only on the 
momentum in absence of an external field. At this stage, k B T is a free param- 
eter and only in the derivation of the equation of state it will be identified 
with the physical temperature. 
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We are able now to evaluate explicitly all other thermodynamic variables and 
provide a complete macroscopic description of a relativistic system at the 
equilibrium. Let us first calculate the probability density defined as 

n = N"U lt = ± f^pfUpfM. (5) 



Since n is a scalar, it can be evaluated in the rest frame where £/ M = (1, 0, 0, 0). 
Setting r = p°/k B T and z = m/k B T, the above integral can be written as 



n=|(A: B T) 3 Jdr(T* re? . (6) 



Let us introduce the g-modified Bessel function of the second kind as follows 

oo 



then, by means of a partial integration of Eq.(6), the particle density can be 
cast into the compact form 

n=—m 2 k B TK 2 (q,z). (8) 

Z 9 



Considering the decomposition of the energy- momentum tensor [9]: T^ v = 
eU^U" — pA^, where e is the energy density, p the pressure and A M1/ = 
_ ^ the equilibrium pressure can be calculated as 

P = -\t^A, u = ~J ^P^Mp) , (9) 



and can be expressed as 

p = -m*(k B T) 2 K 2 (q,z). (10) 

Z 9 



Comparing Eq.(8) with Eq.(10), we obtain 

p = nk B T, (11) 



which is the equation of state of a perfect gas if we identify T as the physical 
temperature of the system. 
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We proceed now to calculate the energy density e as 

1 f d 3 p 



e = T^Uu = - I -J(p»U»ff e \{p) . (12) 



Inserting the previously defined variables r and z and using the definition in 
Eq.(7), we obtain 



4vr 4 
e = — m 



3 K 2 (q,z) Ki(q,z) 

z 2 z 



(13) 



Thus the energy per particle e = e/n is 



e = 3 *- T+m tj^' (14 » 



which has the same structure of the relativistic expression obtained in the 
framework of the equilibrium Boltzmann-Gibbs statistics [9]. 

For a system of particles in degenerate regime the above classical distribution 
function (4) has to be modified by including the fermion and boson quan- 
tum statistical prescriptions. For a dilute gas of particles and/or for small 
deviations from the standard extensive statistics (q ~ 1) the equilibrium dis- 
tribution function, in the grand canonical ensemble, can be written as [10] 

n(k, /i) = (1 - ^ _ _ ^p^ri) ± i , (15) 



where the sign + stands for fermion and — for boson particle. All the previous 
results can be easily extended to the case of quantum statistical mechanics. 



3 Nonextensive statistics in hadronic matter and QGP 



The motivation of the importance of non-standard statistical effects in nu- 
clear and high energy physics lies in the fact the extreme conditions of density 
and temperature in high energy nuclear collisions give rise to memory ef- 
fects and long-range color interactions. These conditions imply the presence 
of non-Markovian processes in the kinetic equation affecting the thermal- 
ization process toward equilibrium as well as the standard equilibrium dis- 
tribution [11,12]. A rigorous determination of the conditions that produce a 
nonextensive behavior, due to memory effects and/or long-range interactions, 
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should be based on microscopic calculations relative to the parton plasma 
originated during the high energy collisions. At this stage we limit ourselves 
to consider the problem from a qualitative point of view. However, it is note- 
worthy to notice that in proximity of the hadronic-QGP phase transition, 
non-perturbative QCD calculations become important. Only a small number 
of partons is present in the Debye sphere: the ordinary mean field approxima- 
tion of the plasma is no longer correct and memory effects are not negligible. In 
addition, we observe that at high density the color magnetic field remains un- 
screened (in leading order) and long-range color magnetic interaction should 
be present at all temperatures. 

From the above considerations it appears reasonable that in regime of high 
density and temperature both hadronic and quark-gluon EOSs must be af- 
fected by nonextensive statistical effects. In the next two subsections, we will 
study the two EOSs separately on the basis on the previously obtained rela- 
tivists thermodynamic relations. 



3. 1 Nonextensive hadronic equation of state 



Concerning the hadronic phase we use a relativistic non-linear model based 
on the interacting many-particle system consisting of nucleons, isoscalar and 
isovector mesons (a, uo and p)[13]. On the basis of the Eqs.(2), (9) and (12), 
the pressure and the energy density can be written as 



£2 



1 3 1 L 4 i 1 22,1 22 

-aa - -bo- + -m^uj + -m p p , 

d 3 k _ I 



(16) 



t=n,p 



—E;(k)[n«(k,ri)+n«(k,-ri)} + -mla 2 



+ I a(T 3 + 1^4 + l m 2 ^2 + l m 2 p2? 



(17) 



where n(k, fii) and n(k, —fit) are the fermion particle and antiparticle distribu- 
tion (15). The nucleon effective energy is defined as E* = ^k 2 + M* 2 , where 
M* = M N — g a a. The effective chemical potentials fx* are given in terms of 
the vector meson mean fields fii — fi* — g^u =)= g p p (— proton, + neutron), 
where pi are the thermodynamical chemical potentials fii = de/dpi. At zero 

temperature they reduce to the Fermi energies Ep,- L = \J k 2 Fi + M* 2 and the 
nonextensive statistical effects disappear. The isoscalar and isovector meson 
fields (a, uj and p) are obtained as a solution of the field equations in mean 
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field approximation and the related couplings (g a , and g p ) are the free 
parameter of the model [13]. Finally, The baryon densities ps are given by 

/d 3 k 
j^[n(k,ri)-n(k,-ti)\, (18) 

where 7 is the spin/isospin multiplicity. 

In Fig. 1, we report the obtained hadronic EOS: pressure as a function of 
the baryon number density (in units of the nuclear saturation density p = 
0.16 fm~ 3 ) for different values of q. Because of previous phenomenological 
studies in heavy-ion collisions have brought to values of q greater than unity 
[2], we concentrate our analysis to q > 1. The results are plotted at the 
temperature T = 100 MeV and at fixed value of Z/A = 0.4. The range of the 
considered baryon density and the chosen values of the parameters correspond 
to a physical situation that can realized in the recently proposed high energy 
heavy- ion collisions experiment (like in Au — Au collisions) (see, for example, 
Ref.[15]). 
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Fig. 1. Hadronic equation of state: pressure versus baryon number density (in 
units of the nuclear saturation density po) for different values of q. 

3.2 Nonextensive QGP equation of state 

In the simple model of free quarks in a bag, the pressure, energy density and 
baryon number density for a Fermi gas of quarks in the framework of nonex- 
tensive statistics (see Eqs.(l), (2), (9) and (12)) can be written, respectively, 

as 

P = E l^Jk^[n<(k,p f )+n<(k,-H)]k?dk-B, (19) 

f=u,d o 
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oo 

e = E ^ / Vf) + -H)\k 2 dk + B, (20) 

f=u,d o 

oo 

P = E 3 ^ / K*, - -/x/)]A; 2 ^, (21) 

f=u,d o 

where e/ = (k 2 +mj) 1 ^ 2 , n(k, /i/) and n(k,fj,f) are the particle and antiparticle 
quark distributions. The quark degeneracy for each flavor is 7/ = 6. Similar 
expressions for the pressure and the energy density can be written for the 
gluons as massless Bose gas with zero chemical potential and degeneracy factor 
7 9 = 16. Because of the fermion-boson nonextensive distribution (15), the 
results are not analytical, also in the massless quark approximation. Therefore, 
a numerical evaluations of the integral must be performed. Note that a similar 
calculation, only for the quark-gluon phase, was also performed in Ref . [14] by 
studying the phase transition diagram. 

In Fig. 2, we report the EOS for massless quark u, d and gluons, for different 
values of q. As in Fig. 1, the results are plotted at the temperature T = 100 
MeV and at fixed value of Z/A — 0.4; the bag parameter is 5 1//4 =170 MeV. 
In both the figures we can observe remarkable effects in the shape of the EOS 
for small deviations from the standard statistics. 




p/Po 



Fig. 2. Same of Fig. 1 for the case of the quark-gluon equation of state with 



J B 1 1/4 =170 MeV. 



4 Conclusion 



In this work we have studied the relativistic thermodynamic relations and de- 
rived the EOS of a gas of free particle at the equilibrium in the framework of 
the nonextensive Tsallis thermostatistics. The results are applied to obtain a 
consistent generalization of the EOS of strongly interacting hadronic matter 
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and of deconfined QGP. The range of density and temperature chosen are 
physical values estimated in the recently proposed high energy heavy-ion col- 
lisions experiments at finite baryon chemical potential [15]. We find that small 
deviations from the Boltzmann-Gibbs statistics implies a sensible modification 
of the two considered EOSs. A complete discussion on the nuclear physics con- 
sequence lies out the scope of this paper, however we want to outline that such 
a modification of the EOSs can strongly affect different physical properties, 
like, for example, the critical phase transition density, the symmetry energy, 
the nuclear compressibility, connected to experimental observables. 
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